Quantum phases and phase transitions of interacting bosons in a Kagome lattice 
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By using the dual vortex method and analyzing the saddle point structures of the dual gauge 
fields self-consistently in the insulating side, we study quantum phases, especially supersolid phases 
and phase transitions in an extended boson Hubbard model at and slightly away from 1/3 ( 2/3 ) 
filling on a kagome lattice. We find phases and phase transitions in Kagome lattice are dramatically 
different from those in bipartite and triangular lattices in both Ising and Easy-plane limit. Slightly 
away from / = 1/3(2/3), there is no supersolid in the Ising limit, but there are always supersolids 
in the Easy-plane limit where we find a stripe solid and a new type of supersolid which has charge 
density wave- valence bond-superfluid order and also analyze their corresponding stability conditions. 
Implications on QMC simulations with both nearest neighbor and next nearest neighbor interactions 
are given. 



1. Introduction. The Extended Boson Hubbard 
model (EBHM) with various kinds of interactions, at 
various kinds lattices ( bi-partisan or frustrated ) at 
various kinds of filling factors ( commensurate or in- 
commensurate ) is described by the following Hamilto- 
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where n, — h[hi is the boson density, t is the hopping 
amplitude, U,Vi,V2 are onsite, nearest neighbor (nn) 
and next nearest neighbor (nnn) interactions between 
the bosons. The • • • may include further neighbor in- 
teractions and possible ring-exchange interactions. A su- 
persolid in EqnlT]is defined as to have both ofi-diagonal 
long range order < bi >^ and diagonal charge density 
wave in the boson density n^. Possible supersolid states 
in frustrated lattices such as triangular and Kagome lat- 
tices were first investigated in [l|] by spin wave expansion. 
A Kagome lattice ( Fig. la ) has 3 sublattices a,b and c, 
it contains both hexagon and triangle, its dual lattice is 
the dice lattice ( Fig. 2b). From spin wave expansion, the 
authors in [l| found that a SS state is more robust in a 
triangular lattice with only t and Vi terms in Eqn[l]and a 
SS state with \/3x \/3 pattern is stable even at half filling 
/ = 1/2. Indeed, this discovery was confirmed by several 
recent QMC simulations [3|. However, due to too strong 
quantum fluctuations, the authors in [l| found that the 
spin wave expansion does not work anymore in Kagome 
lattice. So a completely different analytical method is 
needed to study possible supersolid phases in a Kagome 
lattice. 

Recently, the EBHM in square lattice at generic com- 
mensurate filling factors / = p/q { p,q are relative 
prime numbers ) were systematically studied in [3| by a 
dual vortex method ( DVM ). However, it was explicitly 
pointed out in [5| that the DVM developed in pi| holds 



only in the featureless superfluid ( SF ) and the valence 
bond solid ( VBS ) side where the saddle point of the 
dual gauge field can be taken as uniform, but it fails in 
the charge density wave ( CDW ) side where the saddle 
point of the dual gauge field can not be taken as uni- 
form anymore. So the DVM developed in 3l| is not self- 
consistent in the CDW side, the CDW order parameter 
constructed in [3| to characterize the symmetry breaking 
patterns in the CDW side are also not self-consistent. So 
special care is needed to choose a correct saddle point of 
the dual gauge field in the CDW side, so a different effec- 
tive action is needed in the CDW side to make the theory 
self-consistent [5|. In [5[, (1) by extending the DVM ex- 
plicitly to the lattice symmetry breaking side by choos- 
ing the corresponding self-consistent saddle points of the 
dual gauge field (2) by pushing the DVM to slightly away 
from commensurate filling factors, the author mapped 
out the global phase diagram of the EBHM on bipartite 
lattices such as honeycome and square lattice near half 
filling and investigate superfluid, solid, especially super- 
solids and quantum phase transitions in the global phase 
diagram in a unified scheme. 

In this paper, I will use the DVM developed in [3|, jj, |5| 
to study quantum phases and phase transitions in the 
kagome lattice at and near filling factor / = 1/3 (Fig. 2 ). 
In a frustrated lattice, the sign of the hopping t in Eqn[T] 
makes crucial differences. It is also well known that the 
physics in Kagome lattice is much more complicated than 
that in a triangular lattice as indicated by the spin wave 
calculations in [l!| . We will also compare our results with 
some previous Quantum Monte-Carlo (QMC) results in 
the Kagome lattice [8|, |9| and also give important im- 
plications to possible future QMC simulations on Eqn[T] 
with both hard core and soft core, with both the nn Vi 
and the nnn V2 interactions. 

2. The effective action. Because of the P-H symmetry 
at [/ = CX3, the results are equally valid near 2/3 filling. 
As shown in y at g = 2, because the lowest dual vortex 
band in the dice lattice is completely flat, so the dual 
vortices are completely localized, so it has to be a super- 



fluid, in sharp contrast to g = 2 on triangular lattice [a|. 
For q = 3, from Table 1 of |4|, we can see there are two 
minima at (0,0) and {—2tt/3,2tt/3). Let's label the two 
eigenmodes at the two minima as ipij = 0,1. The effec- 
tive action invariant under all the MSG transformations 
upto sixth order terms was written down in [Sj. Again, 
inside the superfluid phase, moving slightly away from the 
1/3 filling / = 1/3 corresponds to adding a small mean 
dual magnetic field Sf = f —1/3 in the action derived in 
[8|. Upto sixth order, the action is Csf = Cq + Ci + C2' 
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where A^ is a non-compact U{1) gauge field. 

In fact, the form of this action is exactly the same as 
that in honeycomb lattice at g = 2 first derived in [5|. 
However, the physical meanings of the two eigenmodes 
tfji,l = 0,1 are completely different in the two lattices, 
this fact leads to completely different physics in the two 
lattices. In the following, we focus on the insulating side 
where r < 0. 

3. Ising limit. If v > 0, the system is in the Ising 
limit, only one of the 2 vortex fields condenses. Let us 
take {tpo) 7^ 0, (-^i) = 0. It was found that there is no 
CDW order in [8], but there is a VBS order around ei- 
ther up or down triangles (Fig.l). This is just opposite to 
bipartite and triangular lattices discussed in [3, H, la, Q 
where the Ising limit corresponds to the CDW case. We 
interpret "0;, ^ = 0, 1 as the kinetic energies around the up 



and down triangles Kijk 
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where i,j,k are the three vertices of a triangle, namely, 
themselves are valence bond order parameters, while in 
bipartite and triangular lattices, they are CDW order pa- 
rameters. This interpretation is indeed consistent with 
the phases to be identified in the Easy-plane limit in the 
following. Because the saddle point structure stays the 
same across the SF to the triangle- VBS ( TVB ) tran- 
sition at / = 1/3, EqnE] still holds in the TVB side, 
so the transition is a weak first order one which breaks 
both the U{1) and the Z2 exchange symmetry between 
■00 and 01 [12||. In contrast, in bipartite and triangular 
lattices, in the Ising limit, different saddle points for the 
dual gauge fields need to be chosen [^, [y] . These crucial 
differences are responsible for the absence of supersolid 
in the Ising limit in a Kagome lattice as shown in the 
Fig.l. Slightly away from 1/3 filling, there must be a 
direct frist order transition from the TVB to the super- 
fiuid, there is no TVB supersolid intervening between the 
TVB and the SF in the Fig.l, in sharp contrast to the 
cases in bipartite and triangular lattices. 

4. Easy-plane limit. If u < 0, the system is in the 




tvb\ 


SF 


<V„>^o\ 


<%>=() 


<V,>=o 


\ <^>=() 



(a) 



r<0, v>0 

(b) 



FIG. 1: ( a ) The triangular valence bond ( TVB ) in a 
Kagome lattice at / = 1/3. The thick triangle means one 
boson is moving around the up triangle ( b ) The phase di- 
agram driven by the chemical potential fi. There is a direct 
first order transition from the TVB to the superfluid. There 
is no supersolid intervening between the TVB and the SF. 
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FIG. 2: ( a ) A Kagome lattice has 3 sublattices labeled as 
a, h, c. The CDW-I-VBS in a Kagome lattice at / = 2/3 iden- 
tified in [9|]. The dot means one ( or 1 — a ) boson occupying 
the site. The O inside the hexagon means 3 ( or 3(1 -I- a) 
bosons hopping around the hexagon. ( b ) The saddle point 
of the dual gauge field in a magnetic fleld unit cell in the dice 
lattice. 



easy-plane limit, Eqn[2] is similar to the action in Bi- 
layer quantum Hall systems _11|, the two vortex fields 
have equal magnitude 0/ — |0|e*^' and condense. The 
relative phases can be determined by the sign of ^[5]. In 
this case, the system is expected to have either CDW or 
both CDW and VBS around a hexagon. As explained in 
the introduction, if there is any CDW component, Eqn[2] 
does not work anymore, so can not be used to character- 
ize the symmetry breaking pattern in the insulating side. 
However, the results in [q| can still be used as a useful 
guide to see some signatures of the symmetry breaking 
patterns. In the following, we discuss both cases sepa- 
rately. 

(a) Both CDW and VBS order: The analysis in d 
suggests that this is the case when w > 0, so 9^ = 
6»o - 6*1 = (2n -I- l)7r/3 in EqnO Eqn[T]in Kagome lattice 
with [/ = 00, Vi > was studied by recent QMC Simula- 



tions [9|, |lO| . It was found that the sohd state at exactly 
1/3 fining has both CDW and VB order. The VB order 
corresponds to boson hopping around a hexagon ( shown 
in Fig. 2a ). This is a unique feature of the Kagome lat- 
tices. This feature can be intuitively understood with- 
out any calculation: with very large Vi interaction at 
1/3 filling, if a boson taking sublattice a, then all its 
4 nn which belong to sublattices b and c are excluded, 
so bosons will go to the 4 nnn which still belong to b 
and c instead of a ( Fig. 2a). This is in sharp contrast 
to triangular lattice at 1/3 fiUing in [6|. So bosons may 
move around the hexagon to form some VB order in ad- 
dition to the CDW order ( Fig. 2a). Indeed, the up and 
down triangles in the Fig. 2a, each containing one CDW 
site and one VB, are equivalent. The saddle point for 
the dual gauge fields is < V x Acdw >= 1 — a, < 
V X AvBS >= (1 + a)/2 as shown in Fig.2b. The 
lowest energy dual vortex band in such a saddle point 
is found to be: E{ki,k2) = — ■\/6 + 2A(fci, fe) where 
A(fci, fe) = cosfci— cos(fci — a7r) + cos/i;2 — cos(fc2-l-a7r)-|- 
cos(/ci — fc2) — cos(fci — fc2 — 2a7r), — tt < fci, k2 < tt. It is 
easy to see that there is a symmetry E{ki,k2) — E{k2 + 
7r(l -I- a), ki + 7r(l - a)). For < a < 0.404, there are 
two minima at fci = a7r/2 — arcsini ^^'^,' s ),fc9 = —fci 

-■- / ^2 cos(Q7r) /' ^ -L 

and fci — om 12 + arcsin C^,' \ \ + tt. k2 = —fci . The 
two minima are related by the symmetry. 

When a = 1/3 < 0.404, it reduces to q — 3 case where 
there are two minima at (0,0),(— ^,^) which is the 
two minima at the superfluid side in Eqn[21 For < 
a < l/3~, let's still label the two minima by tpi, I = 0,1. 
The a — > l/3~ case corresponds to approaching to the 
transition from the CDW-I-VB to the SF in Fig. 3a from 
the left hand side. Then the MSG inside the SF state 
Ti,T2,RTr/3,Ii is reduced to Ti, T2, i?^, /s^+Sa inside the 
CDW-HVBS state in Fig.2b. The effective action slightly 
away from 1/3 invariant under the MSG transformations 
inside the CDW-fVBS state upto sixth order terms is 
^CDW+VBS ^ Co + Ci + £2- 

Co = 5]|(a^-z<^^)^,p+r>d' 
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where AY_^^ is a non-compact C/(l) gauge field. Note 
that A'r^^^ is always massive, so was already integrated 
out in Eqn[31 So the density fluctuations on the dot sites 
in Fig. 2a are suppressed. Because the system is already 
in the CDW-I-VBS, so we expect f < and v <0. Then 
we can set tpi = |V'|e*^'. The relative phases can be de- 
termined by the sign of w > which has the same sign 
as w; > 0. It is easy to observe that due to the same 
structure of the saddle points on the SF side Eqnl^] and 



on the CDW-VBS side EqnlH the SF to the CDW-HVB 
transition driven by r along the horizontal axis in Fig. 3a 
is very weakly first order. This is in sharp contrast to 
the strong 1st order SF to CDW transition in a triangu- 
lar lattice at / = 1/3 studied in |6|]. It is also different 
from the SF to the VBS transition at f — 1/2 in honey- 
comb and square lattices where the same action can be 
used on both side. Although the QMC in Ref. ^IQ] found 
the CDW-VB solid to the SF at / = 1/3 is a 2nd order 
transition through a possible dcconfincd quantum criti- 
cal point, the QMC in Ref.[£^] found a weakly 1st order 
transition from double-peaked histogram of the boson ki- 
netic energy. Our results from the self-consistent DVM 
explains the underlying mechanism of the weak first order 
transition. 

When moving from 1/3 to 1/2 where the boson has to 
be a SF, Eqnl3]can be rewritten as: 



C = i^d.e^ -Al^^r + ^{e,.,d.Ar'- 2nd fS,^)' 
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where d± = 6*0 ± ^1 ■ 

The SS slightly away from 1/3 filling has the same 
CDW and VB order as the solid at 1/3, so let's call 
this kind of novel SS as CDW-VB-SS. Then the transi- 
tion from the CDW-t-VB to the CDW-VB-SS transition 
driven by the chemical potential in EqnU] can be sim- 
ilarly discussed as the VBS to the VB-SS in bipartite 
lattices in [5|, la, lit , so it is is the same universality class 
as that from the Mott insulator to the SF. The CDW- 
VB-SS to the SF transition is 1st order. Although the 
CDW-VB-SS may not be stable in the hard core case, it 
should be stable in the soft core case. Of course, there 
is no P-H symmetry anymore in the soft core case. This 
study may inspire more accurate QMC to search for this 
novel CDV- VB-SS state, especially in the soft-core case. 

(b) Only CDW order: We expect that this is the case 
when w < 0, so 0^ ^ 00 - 01 ^ 2n7r/3 in EqnH If 
we add very strong nnn interaction V2, the bosons will 
go to the nnnn which belong to a again ( Fig. 2a), then 
the bosons will simply take one of the 3 sublattices such 
as a to form a CDW alone ( Fig. 3b). Indeed, the up 
and down triangles in the Fig. 2a, each containing one 
occupied site and two empty sites, are equivalent. In 
the stripe CDW state in Fig. 3b, a different saddle point 
where < V xA^ >= 1 — 2a for the sublattice a and < V x 
A*"^ >= a for the two sublattices b and c should be used. 
The a ^ etc < 1/3 limit corresponds to approaching the 
stripe CDW to the SF transition from the stripe CDW 
side in Fig. 3b, while the a — *■ 0"*" limit corresponds to 
t/Vi -^ limit in Fig. 3b. It is easy to see that there is 
only one vortex minimum ^tc in such a staggered dual 
magnetic field with a < ac < 1/3, so the effective action 
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FIG. 3: ( a ) Zero temperature phase diagram in Kagome 
lattice U = oo, Vi > 0. In the soft core case U < oo, 
there maybe a narrow window of CDW-VB-SS sandwiched 
between the CDW-VB sohd and the SF. In the hard core 
case studied in [y, [l0|, the CDW-VB-SS is unstable, there 
is a direct weakly first order transition from the CDW-VB 
solid to the SF. (b) Very similar phase diagram also works for 
U = CO, Vi > 0, V^2 > 0. In this case, Vi in (a) need to be 
replaced by Vi -|- V2, the solid at 1/3(2/3) is simply a CDW 
where the bosons occupying ( un-occupying ) one of the 3 sub- 
lattices a, b, c in Fig. 2a. The SS is simply a stripe SS where 
the superfluid stiffness ptc ^ pba- The Stripe-SS will surely 
be stable even in the hard core case. As shown in the text, 
the first order transition in (a) is a weakly first order one, but 
strongly first order one in (b). Thin (thick) line is 2nd (first 
) order transition. 



inside the stripe CDW state is: 
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where the vortices in the phase winding of ipbc should be 
interpreted as the the boson number. Because of sharp 
change of the saddle point from Eqnl2] in the SF side to 
Eqnl5]in the stripe CDW side, the transition from the SF 
to the stripe CDW along the horizontal axis in Fig. 3b is 
likely to be a strong first order. 

Slightly away from 1/3, there is a stripe SS as shown in 
Fig. 3b. In the direct lattice Fig. 2a, even if taking away 
the sublattice a, the bosons can still move easily along the 
chain bcbc..., so they need only overcome the barrier along 
the chain baba... to achieve an effective hopping along 
baba.... So the superfluid stiffness along the be chain is 
larger than that along the ba chain. Because the nn 
and the nnn neighbors are similar as can be seen from 
Fig. 2a, the vertical axis in Fig. 3b need to be replaced 
by m/(^i + Vi), the horizontal axis stay as t/Vi at fixed 
large V2 ( Fig. 3b ). While the conclusion achieved in 
[J| that the system has to be a superfluid state at 1/2 
remains robust. So we conclude that when moving from 
1/3 to 1/2 driven by the chemical potential along the 
vertical axis in Fig. 3b, the system goes CDW to stripe- 
SS to the SF transition. The first is the second order, the 
second is a first order transition ( Fig. 3b). Again, longer- 



range interactions favor the description in terms of the 
dual vortices and the stability of SS, so the Stripe-SS 
will surely be stable even in the hard core case ( Fig. 3b 
). We suggest that it is important to do QMC with U — 
00, Vl > 0, ^2 > to test Fig.3b. 

5. Conclusions. By using the DVM developed in 
[3|, |5|, |6[, we studied quantum phases and phase tran- 
sitions in an extended boson Hubbard model at and 
slightly away from 1/3 ( 2/3 ) filling on a kagome lat- 
tice. We find that in the insulating side, at the com- 
mensurate fillings such as 1/3, one has to choose the 
corresponding saddle point structure of the dual gauge 
fields self-consistently to construct effective actions on 
this side. We find quantum phases and phase transi- 
tions in the Kagome lattice in both Ising and easy plane 
limit are dramatically different from those in bipartite 
and triangular lattice studied previously [2, IS, la, la, [fl] . 
We interpret the two vortex order parameters tjji,l = 0,1 
as the kinetic energies around the up and down triangles 
in a Kagome lattice. In the Ising limit, the SF to the 
triangle valence bond order (TVB) is a weak first order 
one. For the first time, the Ising limit of the same action 
can be used to describe a quantum phase transition in a 
concrete model. In the easy-plane limit, we find that the 
CDW-I-VB to the SF transition at commensurate filling 
/ = 1/3 in Kagome lattice is a weakly first order one. 
This finding resolves some outstanding puzzles observed 
in the QMC in a kagome lattice [9|, llO| . When pushing 
the effective action slightly away from the commensurate 
fillings, we find that in the Ising limit, there can only 
be a direct first-order transition from the triangular va- 
lence bond (TVB) to the SF, so there is no immediate 
TVB supersolid intervening between the TVB and the 
SF. However, in the Easy-plane limit, different kinds of 
supersolids are generic stable states slightly away from 
1/3 ( 2/3 ) fiUing. In addition to a stripe supersolid ( 
stripe-SS) , we also find a new kind of supersolid: CDW- 
VB supersolid. The stripe-SS should be stable even in 
the hard core limit, while the CDW-VB-SS maybe stable 
only in the soft core limit. We also made implications 
on the available and future QMC simulations on EBHM 
with both nn and nnn interactions in a kagome lattices. 

I am indebted to Longhua Jiang for technical as- 
sistance in finding the saddle point structures in the 
Kagome lattice. 
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